Qubit-oscillator systems in the ultrastrong-coupling regime and their potential for 

preparing nonclassical states 
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We consider a system composed of a two-level system (i.e. a qubit) and a harmonic oscillator 
in the ultrastrong-coupling regime, where the coupling strength is comparable to the qubit and 
oscillator energy scales. Special emphasis is placed on the possibility of preparing nonclassical 
states in this system. These nonclassical states include squeezed states, Schrodinger-cat states and 
entangled states. We start by comparing the predictions of a number of analytical methods that 
can be used to describe the system under different assumptions, thus analyzing the properties of the 
system in various parameter regimes. We then examine the ground state of the system and analyze 
its nonclassical properties. We finally discuss some questions related to the possible experimental 
observation of the nonclassical states and the effect of decoherence. 



I. INTRODUCTION 

The two-level system (or qubit) and the harmonic os- 
cillator are the two most basic, and perhaps most often 
studied, components of physical systems. The paradigm 
of a qubit coupled to a harmonic oscillator has also 
been analyzed by numerous authors over the past few 
decades [l|, Q • Physical systems that can be described 
by this model include natural atoms coupled to optical 
or microwave cavities [3] , superconducting qubits coupled 
to superconducting resonators [^-01, quantum dots or 
Cooper-pair boxes coupled to nanomechanical resonators 
electrons interacting with phonons in a solid fol] and 
some models of chaotic systems 10]. 

In the early work on cavity quantum electrodynam- 
ics (QED) in atomic systems, the achievable atom-cavity 
coupling strengths were smaller than the atomic and cav- 
ity decay rates, usually limiting observations to only indi- 
rect signatures of the theoretically predicted phenomena. 
Recently, the strong-coupling regime, where the coupling 
strength is larger than the decay rates in the system, has 
been achieved [Tl|. In addition to atomic systems, the 
strong-coupling regime has been achieved in supercon- 
ducting circuit-QED systems and superconducting- 
qubit-nanomechanical-resonator systems are approach- 
ing this regime 0- In fact, superconducting systems 
are suited for achieving the so-called ultrastrong-coupling 
regime, where the qubit-oscillator coupling strength is 
comparable to the qubit and oscillator energy scales [l3] ■ 
One can expect to find new phenomena in this regime 
that are not present in the weak or moderately strong 
coupling regimes. Indeed there have been a number of 
theoretical studies on this system analy zing some of its 
rich static and dynamical properties |13l4l8l|. 

One reason why superconducting systems are well 
suited for the implementation of qubit-oscillator experi- 
ments is the flexibility they allow in terms of designing 
the different system parameters. For example, in the 
two earliest experiments on circuit QED, Chiorescu et 



al. 3 1 used a low- frequency oscillator, while Wallraff et 
al. [J| realized a resonant qubit-oscillator system. Sub- 
gigahertz qubits have also been realized in recent exper- 
iments [19|, and there should be no difficulty in fabri- 
cating high-frequency oscillators. Therefore, all possible 
combinations of qubit and oscillator frequencies are ac- 
cessible, in principle. One advantage of superconducting 
qubits over natural atoms is the additional control as- 
sociated with the tunability of essentially all the qubit 
parameters (20j . as will be discussed in more detail be- 
low. This tunability contrasts with the situation encoun- 
tered with natural atoms, where the atomic parameters 
are essentially fixed by nature. This advantage can be 
seen clearly in the recent experiments where Fock states 
and arbitrary oscillator states were grejDared in a super- 
conducting qubit-oscillator system 121, 22^. We shall see, 
however, that the additional controllability comes at the 
price of having to deal with additional coupling chan- 
nels to the environment, and this unwanted coupling can 
increase the fragility of nonclassical states. 

In this paper we present analytical arguments and nu- 
merical calculations pertaining to the strongly coupled 
qubit-oscillator system from the point of view of the po- 
tential for preparing nonclassical states in this setup. 
These states include squeezed states or superpositions 
of macroscopically distinct states (i.e. Schrodinger-cat 
state) in the oscillator, as well as qubit-oscillator entan- 
gled states Q. In this study, we shall consider all the 
different combinations of qubit and oscillator frequencies. 
We shall also analyze in some detail the effect of the tun- 
ability in the qubit parameters on the behaviour of the 
system. 

The paper is organized as follows: In Sec. |lT]we intro- 
duce the Hamiltonian that we shall use throughout the 
paper. In Sec. IIIII we discuss various analytical meth- 
ods that can be used to study the system under different 
assumptions, and we compare the predictions of these 
methods. In Sec. IIVI we present results of numerical 
calculations that demonstrate the properties of the en- 
ergy eigenstates of the system, including the nonclassical 
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properties of the ground state. In Sec. |V]we discuss the 
possibihty of preparing and detecting the three types of 
nonclassical states of interest. In Sec. IVII we discuss the 
effect of decoherence on the robustness of nonclassical 
states. Section IVlIl contains some concluding remarks. 



II. HAMILTONIAN 




qubit 



coupling oscillator 



FIG. 1: Schematic diagram of the system under consideration. 
A qubit with energy separation _Eq between its ground and 
excited states {\g) and |e)) is coupled to a harmonic oscillator 
with characteristic frequency loo. The coupling strength is 
given by g or A, depending on the language used to describe 
the oscillator. 

The system that we consider here is a qubit coupled 
to a harmonic oscillator, as illustrated in Fig. [H Rather 
than worry about deriving the model from a microscopic 
description of an electric circuit (see e.g. Refs. [H, [23; 
[25}). we shall assume that the description of the system 
as being composed of these two physical components with 
a coupling term of the standard form is an accurate de- 
scription of the system. 

The Hamiltonian of the system is given by: 



H — Ha + Hho + Hw 



(1) 



where 



^ho = 
Hint 
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gxaz 



(2) 



and tr^ are the usual Pauli matrices (with CTz It) === 
It) i^z It) = — |t))j and X and p are the position and 
momentum operators of the harmonic oscillator. The 
parameters A and e are the gap and bias which charac- 
terize the qubit, m is the oscillator's effective mass, ujq 
is the oscillator's characteristic frequency, and g is the 
qubit-oscillator coupling strength. Note that in contrast 
to atomic cavity QED systems, where e = 0, this param- 
eter is easily tunable in present-day experiments using 
superconducting qubits. We shall therefore treat e as 
a tunable parameter (It is worth noting here that most 
past studies on the ultrastrong-coupling regime have fo- 
cused on the case e = 0; however, see Ref. [13) • For 
definiteness, we shall take A and g to be positive. 



It is convenient for some calculations to express the 
oscillator Hamiltonian using the creation (a^) and anni- 
hilation (a) operators: 



a' 
X = 

P = 



X + iP 
X -iP 
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(3) 



which give 
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The coupling strength can therefore be quantified either 
through g or X. 

We shall refer to the eigenstates of Hq as the qubit 's 
ground and excited states, denoted by \g) and |e), keep- 
ing in mind the caveat that this identification becomes 
less meaningful for strong qubit-oscillator coupling. The 
energies of these two states are ztEq/2, where Eq = 
VA^ -I- e^. It is also useful to define an angle 9 that char- 
acterizes the relative size of the and terms in the 
qubit Hamiltonian: tan 6' = e/A. The eigenstates of ifho 
are given by |n), where n — 0, 1,2, with energies given 
by nhujQ (up to the ground state energy ?ia;o/2, which we 
ignore from now on). The integer n represents the num- 
ber of excitations, to which we shall refer as photons, in 
the oscillator. 



III. COMPARISON BETWEEN DIFFERENT 
ANALYTICAL METHODS 

In this section we describe some analytical methods 
that can be used to determine the properties and be- 
haviour of the system based on different assumptions 
(which are valid in different parameter regimes), and we 
compare the predictions of the different methods. 



A. Weak coupling 

The simplest limit is probably the weak-coupling limit 
12], where A ^ Eq.fiujQ. Strictly speaking, one also 
needs to consider the number of photons in the oscillator 
when determining whether the weak-coupling condition 
is satisfied. However, since in this paper we focus on 
a system that remains close to its ground state, we as- 
sume a small number of photons in the oscillator. In the 
weak-coupling limit, one can think of the qubit and os- 
cillator as being well-defined, separate physical systems 
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that interact weakly and can exchange excitations with 
one another [ij- 

In the Hmit of smah A, one can treat Hmt as a small 
perturbation in the total Hamiltonian. The energy eigen- 
states in the absence of this perturbation are given by 
\n,g) — \n) (g) \g) and \n,e) = \n) (g) |e), with energies 
nhLjQ ± Eq/2 (recall that we ignore the TiujqI2 term in 
the oscillator's energy). 

When there are no degeneracies in the non-interacting 
system (i.e. in the Hamiltonian given by iJq + i?ho), the 
addition of the perturbation iJjnt has only a small effect 
on the behaviour of the system. This perturbation only 
slightly modifies the energy levels and eigenstates of the 
Hamiltonian. 

When _Eq w ?iwo, the states \n,g) and |n — l,e) are 
nearly degenerate (note that there is one such pair of 
nearly degenerate states for each value of n), and the 
perturbation term couples them. In particular, the rele- 
vant matrix elements are given by (n, g\ H-mt \n — l,e) — 
X^/ncos9, {n,g\Hint\n,g) = {n~l,e\Hint\n-l,e) = 
0^ and (n-l,e|(i?q -1- Hi,^) \7i - 1 , e) - {n,g\{Hq + 
Hho) \n,g) ~ Eq — HajQ. In other words the effective 
Hamiltonian that one needs to consider is given by 



f 5/2 Xy/Ticose 
[\y/Ecose --5/2 



(5) 



where 5 = Eq — Hloq. Using this Hamiltonian, one can 
analyze the behaviour of the system. In particular, when 
Eq = huiQ, an excitation oscillates back and forth between 
the qubit and oscillator with frequency 2A-\/ncos 9, which 
is commonly referred to as the Rabi frequency. 

Degeneracies also occur when Eq = khujQ, with k being 
any integer. In this case, one can go to higher orders in 
perturbation theory and obtain analytic, though some- 
times cumbersome, expressions describing the properties 
and dynamics of the system. We shall not go any further 
in analyzing this situation here[23|. 

Note that the same results as those explained above 
(for the case Eq « hioa) can be obtained by taking the 
term Hint in the Hamiltonian and replacing it by its 
rotating-wave-approximation (RWA) form: 



-ffint,RWA = A cos 6* (aa+ 



(6) 



where a± are the qubit raising and lowering operators 
((7+ |f;) = |e) etc.). This approximation therefore ignores 
the so-called counter-rotating terms in Hint, which are 
proportional to a^CT-f and aa^, as well as a term propor- 
tional to {a+a''){\e) {e\ — \g) {g\) that appears when e ^ 0. 
These terms would change the number of excitations in 
the system, thus mixing states that have a large energy 
separation (assuming A ^ Eq,hujo), and energy conser- 
vation suppresses such processes. The RWA therefore 
approximates the original Hamiltonian by one where the 
state \n,g) is coupled only to the state \n — l,e), which 
would lead to exactly the same algebra and results men- 
tioned above. Some of the recent studies on ultrastrong 



coupling have analyzed the effects of the counter-rotating 
terms on the system dynamics [l^ [l^ ■ 



B. High-frequency, adiabatically adjusting 
oscillator 



The next limit that we consider is that where the os- 
cillator's characteristic frequency loq is large compared 
to the qubit's energy splitting (i.e. hujQ ^ Eq) and also 
compared to the coupling strength {Hloo ^ A). In this 
case one can say that the oscillator remains in its initial 
energy eigenstate (i.e. ground state, first excited state, 
etc.), and this state follows adiabatically any changes in 
the qubit's state. This case was analyzed theoretically in 
Refs. Bill. 

The procedure for adiabatically eliminating the high- 
frequency oscillator from the problem is straightforward. 
One starts by noting that the qubit is coupled to the 
oscillator through the operator (t^. As a result, one can 
think of the oscillator as always monitoring the qubit 
observable and adjusting to be in the instantaneous 
energy eigenstate that corresponds to that value of ct^ 
(Note here that if the qubit is in a superposition of two 
different az states, each part of the superposition - with 
a well-defined value of az - will have the oscillator in the 
corresponding energy eigenstate). 

We therefore start by assuming that the qubit has a 
well-defined value of Uz, equal to ±1. The oscillator now 
feels the effective Hamiltonian (calculated from /fho and 

-ffint): 



H 



ho. off 



hujoo^a ± X (a + a^) 



(7) 



This Hamiltonian corresponds simply to the original os- 
cillator Hamiltonian with a constant force term applied 
to it. This force term can be eliminated using the trans- 
formation 



a' — a ± 



X 



which gives 



ho.eff 



A^ 

hujoo.''' a' — - — . 



(8) 



(9) 



The above steps can also be carried out in the language 
of the operators x and p: 
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+ —mujQx'^ ± gx 




x' 


= x± 
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-ffho.cff 
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(10) 
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The energy levels of the oscillator are given by nfkuo — 
}?/{hujQ), independently of the qubit's state. There will 
therefore not be a qubit-state-dependent energy that we 
need to take into account when we turn to analyzing the 
behaviour of the (slow) qubit. The oscillator's energy 
eigenstates, however, are slightly dependent on the state 



of the qubit. In particular, 



(11) 



and 



(n<T,=+i|TO<,^=-i) = 



g-2AV(?i"o) 



huja J 



2X 

hujQ 



I n! T va — n 



I rn\_ T n — m 



( 2A 



2\ 



m > n 



m < n, 



(12) 



where Snm is the Kronecker delta, and Ll are the associ- 
ated Laguerre polynomials. 

Having obtained the states of the high-frequency os- 
cillator and the properties of these states, one can now 
turn to the slow part of the system, namely the qubit. We 
take any given value for the index ri, which specifies the 
oscillator's state, and we use it to construct an effective 
qubit Hamiltonian for that value of n. Since there are two 
qubit states for each value of n, the effective Hamiltonian 
will be a 2 X 2 matrix operating in the space defined by 



the states 



I ri,t^, (We use the tildes in order to 

emphasize that the oscillator's state is different from the 
state |n) of the free oscillator). The four relevant matrix 
elements can be calculated straightforwardly as 



Ha 



H, 



A 



Ha 



q 



-2\^ / (huo f j^O 



\fioJoJ 



(13) 



The qubit is therefore described by the effective Hamil- 
tonian 



e A 

A -e 



(14) 



where A is given by Eq. ([T^ and can be thought of as 
the renormalized value of the gap A. 

The exponential and Laguerre-function factors are 
both slightly smaller than one for small values of A/ (huto). 
The qubit therefore experiences a small reduction in the 
coupling (or 'tunnelling') between the states |t) and 
in the weak-coupling limit (A Tt,ujq). This decrease in 
the renormalized value of A can be understood in terms 
of the qubit having to 'pull' the oscillator with it as it 
tunnels between the states |t) and ||), which would slow 
down the tunneling process. Note that the renormalized 
gap depends on the number of photons in the oscillator, 
which can lead to beating dynamics and other interest- 



ing phenomena that occur when several values of n are 
involved jl4 il5i]. 

If we keep increasing A/ (hujo), without worrying about 
satisfying the condition A ?iwo, we find that the La- 
guerre polynomial and therefore the renormalized qubit 
gap vanish at n different A/(?iwo) values (For example, 
for n = 2 there are two values of A for which the renor- 



malized gap vanishes). At these points, the states 



and 71,4^ are completely decoupled. Apart from this 

feature, the renormalized gap decreases as a Gaussian 
function with increasing \/(Tiujq). Note that while in- 
creasing A there is no point that can be seen as a 'critical 
point' with a sudden change in behaviour. This situation 
contrasts with what happens in the two calculations that 
we shall discuss in the next two sections, and the meaning 
of the term critical point will become clearer there. 

One might expect that the adiabatically-adjusting- 
oscillator approximation would break down when A is 
comparable to, or larger than, Tiuq (meaning that hu^ is 
not the largest energy scale in the Hamiltonian) . The ar- 
guments in the previous paragraph might therefore seem 
of little significance. It turns out, however, that the re- 
sults discussed above hold, even when A > fiujQ. The rea- 
son why the approximation of an adiabatically adjusting 
oscillator is still valid in this case is that even though 
large changes in the oscillator's states can occur when 
the qubit's state changes, these large changes involve very 
slow processes that are governed by the renormalized gap. 
The oscillator can therefore adjust adiabatically to these 
slow processes. In other words, the condition that Tiuq 
be the largest energy scale in the Hamiltonian is a suffi- 
cient but not necessary condition for the validity of the 
adiabatically-adjusting-oscillator approximation. 



C. High-frequency, adiabatically adjusting qubit 

We now take the limit where i5q is much larger than 
both TiLOo and A (Some analysis of this case was given in 
Ref. |15,]). Similarly to what was done in Sec. IIII Bl we 
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now say that the qubit remains in the same energy eigen- 
state (ground or excited state), and this state changes 
adiabaticahy fohowing the dynamics of the slow oscilla- 
tor. We therefore start by finding the energy eigenstates 
of the (fast-adjusting) qubit for a given state of the (slow) 
oscillator. Since the interaction between the qubit and 
the oscillator is mediated by the oscillator's position op- 
erator a;, we start the calculation by assuming that x has 
a well-defined value and treat the effective Hamiltonian 
(obtained from iJq and -ffint): 



In the limit Eq ^ g\x\ for the relevant values of x, the 
effective potential in Eq. (jl7p can be approximated by: 



Kft(x) 



where 



2 ° V ^ 

1 ~2 ( ^ £.9 
-main x =F — _„ ^ 

2 ° V mColEq 



egx 
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VA2 



(18) 



q,Gff 
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2 ^ 



2 



gxa^ 



(15) 



(Dn = ± 2- 



The eigenvalues and eigenstates of this Hamiltonian are 
given by: 



\9x) 
E'q,2\x 

|ex) 
tan if 



-^VA2 + (e-2.gx)2 
cos I It) + sin I i;) 

sin I It) -COS I ID 
A 



2gx 



(16) 



We can now take these results and use them to analyze 
the behaviour of the oscillator. We note here that since 
the variable x appears inside the square-root in the above 
expressions, the operators x and p lead to a more trans- 
parent analysis than the operators a and a^. We there- 
fore use the operators x and p for the remainder of this 
calculation. 

Since the qubit's energy depends on the oscillator's po- 
sition X, the oscillator's effective potential now acquires 
a new contribution (which depends on the qubit's state): 



1 



1 



Vcsix) = ^mujlx^ ± + (e 



2gxf 



(17) 



The plus sign corresponds to the qubit being in the ex- 
cited state, and the minus sign corresponds to the qubit 
being in the ground state. In addition to the above effect 
of the qubit on the oscillator, the qubit's state changes 
as the oscillator's position changes, and the oscillator's 
kinetic-energy term will also be modified in principle (this 
effect is similar to the renormalization of A encountered 
in Sec. IIII B|) . However, for a sufficiently high-frequency 
qubit, changes in the qubit's state will be small (see Ap- 
pendix A), and consequently the change in the kinetic- 
energy term can be neglected. 

We now note that the effective potential in Eq. (fT7|) 
is no longer a harmonic potential. The second term de- 
scribes one of the two branches of a hyperbola, depending 
on the qubit's state. It will therefore not be possible to 
derive general analytical results, and we have to start 
considering some special cases. 



mE„ 



(19) 



The oscillator's effective potential is modified in two ways 
depending on the qubit's state. Firstly, the location of 
the minimum is shifted to the left or right by a distance 
proportional to e/Eq = sin 6 (This effect is absent when 
the qubit is biased at the degeneracy point, where e = 0). 
Secondly, the oscillator's frequency is renormalized: Ac- 
cording to Eq. (|19|) . the oscillator's effective frequency is 
increased for the qubit's excited state and reduced for the 
qubit's ground state (This phenomenon is the basis of the 
so-called quantum-capacitance and quantum-inductance 
qubit-readout techniques ^2?]). 

An interesting result appears when one considers the 
case where the qubit is in its ground state and 



2g' 



mujQEq 



> 1, 



4A2 



huJoEc 



> 1. 



(20) 



In this case the renormalized frequency becomes imagi- 
nary. This result signals the presence of a critical point 
above which there is an instability in the system. In par- 
ticular, our expansion of the square-root in Ec^. (1181) is 
no longer valid, the reason being that x would increase 
indefinitely under this approximation and the condition 
-Eq >• g\x\ would be violated. 

The instability obtained above would raise questions 
about the validity of the assumption of an adiabatically 
adjusting qubit above the critical point. Nevertheless, we 
shall not worry about this point now, and we continue 
the calculation. As a first step, we note that Vcf{{x) in 
Eq. PT)) is well behaved at |a;| — > oo. In particular. 



1 



2 2 I 
mUJnX ± 



gx 



when |a;| ^ 



9 9 
(21) 

We can therefore proceed with the calculation using the 
effective potential given in Eq. (fTT)) . 

In order to make a few more statements about the 
case of strong coupling (i.e. above the critical point), it 
is useful to start with the case e = and include a finite 
bias afterwards. When e = 0, the oscillator's effective 
potential takes the form 

(\mul±^)x'±% ,x^A/g 
V,s(x)={ ^ ' (22) 

imwQX^ ± Igxl ,x':^l\/g 
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For the case with the qubit in its excited state (i.e. when 
one has the plus signs in the above expressions), the 
effective potential is a slightly non-harmonic poten- 
tial, and one can expect the oscillator states to look 
more or less like the usual harmonic oscillator states. 
For the case with the qubit in its ground state, and 
when 2g^ / (mujQA) ^ 1, one also has a slightly non- 
harmonic potential. For the qubit 's ground state and 
2g'^ /(mujQA) > 1 (which implies crossing the critical 
point), the oscillator's effective potential is a double- well 
potential. The locations of the minima can be obtained 
by setting dVcs/dx — with Ves given by Eq. ([T7|: the 
minima are located at ±xo, with 



Xo 



A2 



i2,,,4 



(23) 



If one goes well beyond the critical point, the above ex- 
pression reduces to 



Xo 



9 



(24) 



with minimum potential energy [measured relative to 
K=ff(0)] 



2ma;g ' 



and curvature 



dx"^ 



(25) 



(26) 



X — ±Xq 



Note that this curvature is identical to that of the free 
oscillator (i.e. when g = 0). 

One can use the above expressions to estimate the en- 
ergy separation between the ground state and first ex- 
cited state above the critical point. These two states 
will be the symmetric and antisymmetric superpositions 
of the ground states localized around the two minima in 
the double-well potential. The distance between the two 
minima is given by 2xq from Eq. (I24p . and the height 
of the energy barrier separating the two minima is given 
by —Vniin from Eq. ([25|) . Using the Wentzel-Kramers- 
Brillouin (WKB) formula, one finds that the energy sep- 
aration between the two lowest states (and also within 
similar pairs of higher energy levels) is exponential in 
the parameter \/—mVnunXo/^, which is proportional to 
g'^ / {mujQh) , or alternatively / (huJo)'^ . This scaling is 
the same as the one obtained in Sec. IIII Bl 

We now introduce e to the problem. Far below the 
critical point, the effect of e can be obtained easily from 
Eq. (|18|) : the location of the minimum in the single- well 
effective potential is slightly shifted to the left or right. 
More care is required above the critical point, where one 
has the double-well effective potential. In this case, a fi- 
nite value of e breaks the symmetry between the left and 
right wells, thus giving an energetic preference for one of 



the two wells. In order to cause localization in the energy 
eigenstates, e has to be larger than the energy separation 
within one of the energy-level pairs discussed above, i.e. e 
needs to be larger than a quantity that is exponentially 
small in g^ /(mojQh). Clearly, this localization happens at 
smaller values of e as one goes deeper into the bistability 
region. This result means that the superpositions involv- 
ing both wells become increasingly fragile with increasing 
coupling strength. 

Finally we note that above the critical point, one finds 
that the condition Eq ^ g\x\ can no longer be satisfied 
for any of the energy eigenstates. Therefore, one might 
expect that the present approximation cannot be trusted. 
As we discussed in Sec. IIIIBl however, the above results 
hold even when A > Eq. In that case the energy eigen- 
states are either localized close to one of the local minima 
or involve very slow tunneling between the two wells of 
the effective double-well potential. The qubit can adjust 
adiabatically to such slow tunneling processes. It is worth 
mentioning here that when A is the largest energy scale 
in the Hamiltonian, the approximations of this section 
and Sec. IIII Bl are both valid, and either one of the two 
approaches can be used to answer any given question. 



D. Semiclassical calculation 

A semiclassical calculation can go as follows (alterna- 
tive semiclassical calculations can be found in [l^, [3, 
ITsjV The five different variables x, p, ax, Oy and cr^ are 
treated as classical variables whose dynamics obeys the 
Hamiltonian in Eqs. ([ij and ([2]), without the hats. These 
variables obey the constraint C — al. + ay + al — 1. One 
can therefore find the ground state relatively easily by 
minimizing the Hamiltonian under the above constraint. 
Minimizing the function H ~ H — iiC, with /i being a 
Lagrange multiplier, results in the set of equations 



dH 
dx 
dH 
dp 
dH 
dax 
dH 
day 
dH 
da:. 



mujQX + ga^ — 



^=0 



A 

- — - 2^j.ax = 



(27) 



-2^ay ~ 



which are to be solved under the constraint al+al+a'i 



1. The first four equations lead to x 



p = 0, /i = — A/(4cr2:), and ay = 0. The constraint gives 
ax — ±-\/l — a^. One is therefore left with the equation 



± 



Ad, 



2yr 



= 0, 



(28) 
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which can be re-expressed as 



e 

A 



± 



(29) 



The above equation cannot be solved in closed form, in 
general. However, one can make some general state- 
ments about the solution (see Fig. [5]). For the plus sign 
(Fig. [5^), the second term in Eq. ([2^ is a monotonically 
decreasing function that approaches +00 when — — 1 
and approaches —00 when cr^ — >■ 1. There is therefore 
one solution to the equation in that case. It turns out 
that this solution does not correspond to the ground state 
(This fact can be seen as simply a result of comparing the 
energies of the different solutions). The ground state is 
obtained when using the minus sign. In this case, there 
are three possibilities: The first possibility is to have 
< ttioiqA. In this case (Fig. [2t), we find that the 



second term in Eq. (|29p is a monotonically increasing 
function that approaches —00 when ct^ — >■ — 1 and ap- 
proaches +00 when CTz — > 1. There is therefore only one 
solution to Eq. ((29)) . The second and third possibilities 
for the solutions of Eq. (1^^ occur when 2g'^ > toWq A. In 
this case (Figs. [Us and [2k), the second term in Eq. (|^ 
develops a local maximum and a local minimum between 
CTj; = — 1 and tj^ = 1. Depending on the value of e, 
there can be either one or three solutions. In particular, 
when e = 0, the three solutions are given by cr^ = 0, 
which turns out to be an unstable stationary point, and 
az — ±-\/l — (towq A/2g2)2^ which are two degenerate 
ground states (It is easy to verify that this result agrees 
with Eq.l^. 

One can intuitively understand the effect of having a 
finite value of e using the language of Sec. IIII CI For 
e = 0, one has an effective trapping potential for the vari- 
able X, and this potential has the shape of a harmonic- 
oscillator-like single- well potential when 2g^ < towqA 
and a double- well potential when 2g^ > mw^A. This sit- 
uation explains the existence of one ground state when 
2g'^ < mwgA and two degenerate ground states when 
2g^ > towq A. The effect of adding e to the problem is to 
create a tilt in the effective trapping potential; a positive 
value of e favours the negative- a; solution (here we have 
in mind the ground-state solution). If the tilt is weak, 
one has a global minimum in the deeper well and a lo- 
cal minimum in the shallower well. If the tilt exceeds a 
certain critical value, the shallow well is eliminated, and 
one recovers a single-well potential. 

We make a final note on the fact that we started the 
calculation by raising a question related to the ground 
state but found multiple solutions. The reason for this re- 
sult is the fact that Eq. (|27p locates all stationary points, 
and not only the ground state. The calculation there- 
fore identifies both the ground state and also high-energy 
stationary points that are either dynamically unstable 
or dynamically stable but can still relax to lower-energy 
states. 
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FIG. 2: (Color online) Graphical solution of E q. lf29)l (left) 
and the associated effective potentials from Sec. IIII CI (right) 
in the different possible cases. The horizontal lines in the left 
panels represent the first term in Eq. (|29p . and the function 
/((Tz) is the second term in the equation. The circles mark 
the stable solutions of Eq. (|29|) . and the dots in the right 
panels mark the minima in the effective potential of Sec. IIII Cl 
i.e. Eq. (|17[) . The magenta stars in panels (e) and (f) mark an 
unstable stationary point, i.e. a local maximum in the effective 
potential. Panels (a) and (b) correspond to the plus signs in 
Eqs. (|17[) and (|29ll . and panels (c-h) correspond to the minus 
signs. In panels (c) and (d) , the coupling is below the critical 
point, i.e. < mujQEq. In panels (e) and (f), the coupling 
is above the critical point and e is small. In panels (g) and 
(h), the coupling is above the critical point and e is large. 



E. Concluding remarks 

We conclude this section with some remarks on the 
conceptual ideas and predictions of the different analyt- 
ical methods. The standard perturbation-theory proce- 
dure is well suited for the weak-coupling limit. One can 
use it to systematically obtain accurate approximations 
for the energy eigenstates and dynamics of the system. 
The two approximations involving one subsystem, either 
the qubit or the oscillator, adjusting adiabatically to the 
slow dynamics of the other one are based on the con- 
ceptual picture of the separation between different time 
scales in the problem. As formulated above, they in- 
volve only one level of approximation, in contrast to 
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the order-by-order expansion involved in perturbation 
theory. The time-scale-separation-based approximations 
can, however, be constructed formally as the lowest-order 
approximation in a systematic procedure sometimes re- 
ferred to as adiabatic elimination of fast variables [2^ 
or Van Vlcck perturbation theory . The semiclassical 
calculation treats the dynamical variables classically and 
is at first sight not related to any specific approximation 
related to the system parameters. 

The weak-coupling approximation is suited for study- 
ing the excitation-exchange dynamics between the qubit 
and oscillator, but it does not give any hint of an instabil- 
ity in the system. The main result of the adiabatically- 
adjusting-oscillator approximation is the renormalized 
qubit gap. Apart from the oscillatory behaviour in the 
gap, the Gaussian- function decrease at large A/ {hu/o) val- 
ues is a signature of the strong entanglement between the 
qubit and the oscillator in the energy eigenstates. Nev- 
ertheless, no 'critical point', i.e. a point that is associ- 
ated with a sudden change in any of the effective qubit 
parameters (particularly the renormalized gap), is ob- 
tained in that calculation. The adiabatically-adjusting- 
qubit approximation predicts a reduced effective oscilla- 
tor frequency for weak coupling (and assuming that the 
qubit is in its ground state), and a qualitative change in 
behaviour upon crossing the critical point 



hujaEq 

Above the critical point, the energy eigenstates can 
be highly entangled qubit-oscillator states. As in 
the adiabatically-adjusting-oscillator approximation, the 
separation between neighbouring energy levels is found 
to follow a Gaussian-function dependence in the param- 
eter \/(Tiujq). The adiabatically-adjusting-oscillator and 
adiabatically-adjusting-qubit approximations give differ- 
ent predictions regarding the typical value of A at which 
the Gaussian-function decrease in energy separation 
starts: the former gives A ~ Tiloq and the latter gives A ^ 
y/huJoK. The semiclassical calculation produces the same 
critical-point condition as the adiabatically-adjusting- 
qubit approximation. Even though the semiclassical cal- 
culation naturally cannot produce any entangled-state 
solutions, its results can be used as a starting point for 
studying quantum superpositions of the different semi- 
classical solutions. 

One could understand the reason for the absence of a 
critical point in the case of a high-frequency oscillator as 
having to do with the pairing of energy levels. In this 
case, the energy levels form pairs all the way from A = 
to A — )■ cx). In contrast, in the case of a high-frequency 
qubit the low-lying levels are equally spaced for small 
values of A. As A increases, the energy levels start ap- 
proaching each other while remaining equally spaced, a 
situation that corresponds to a decreasing renormalized 
oscillator frequency. At the point where the energy lev- 
els are expected to collapse to a single, highly degenerate 
energy level, they pair up and the different pairs start 



moving away from each other. The energy levels now 
resemble those of an increasingly deep double-well po- 
tential. Thus the energy levels and energy eigenstates 
exhibit two qualitatively different structures below and 
above the critical point. 

It is worth mentioning that the adiabatically- 
adjusting-oscillator and adiabatically-adjusting-qubit ap- 
proximations start with similar, or symmetric, reasoning. 
The asymmetry in the results is mainly due to the dif- 
ferent dependence in the energy levels and energy eigen- 
states of the fast subsystem on the state of the slow sub- 
system. In the case of a fast qubit, the qubit 's energy 
produces the largest effect on the slow oscillator. In the 
case of a fast oscillator, the oscillator's energy does not 
depend on the state of the qubit, and only the changes 
in the energy eigenstates lead to effective changes to the 
behaviour of the slow qubit. 

One might wonder why the results of the semiclassical 
calculation agree with those of the high-frequency-qubit 
approximation. This agreement can be understood by 
noting first that the oscillator has continuous variables, 
such that it is conceivable that certain states will be de- 
scribed to a good approximation using classical variables 
(In this context one can think of coherent states, which 
to a good approximation behave classically). When the 
qubit 's frequency is high, it is also conceivable that the 
qubit 's state, which for example follows the instantaneous 
ground state, can be described by the classical variables 
that specify the instantaneous ground state. In this case, 
one can expect the semiclassical calculation to give good 
results. One can also use this argument to conclude that 
a phase-transition-like singularity will occur in the limit 
hujo/Eq — > 0, where the semiclassical calculation can be 
expected to give exact results. In contrast, in the limit 
of a low-frequency qubit, the dynamics will necessarily 
be described by the coupling of two discrete quantum 
states, and this situation cannot be described well using 
a semiclassical approximation. 

In the context of discussing the phase-transition-like 
bifurcation in this system, it is worth mentioning a re- 
lated s yst em with a true phase transition: the Dicke 
model [301 • If one replaces the single qubit by a large 
number of qubits with equal values of Eq, all coupled 
to the same cavity with the same value of A, then by 
taking the appropriate thermodynamic limit (N — > oo, 
A^iV = A^) one finds a phase transition between states 
similar to those discussed above [sj]. The critical point 
is given by the condition 4A^ = fiuj^Eq, in analogy to the 
critical-point condition discussed in Sees. IIII Gl and UlI Dl 
In contrast to the single-qubit case, however, the phase 
transition now occurs regardless of the relation between 
the qubit and oscillator frequencies. Note that the qubits 
behave collectively as a single large spin in this case (and 
for low-lying states a large spin behaves similarly to a 
harmonic oscillator), such that the entire system can be 
approximated by two coupled oscillators. Note also that 
the semiclassical calculation arises naturally in this case: 
when the effective spin has an infinite number of allowed 
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states, it is natural to make a classical approximation 
where fluctuations in the spin are small compared to the 
total size of the available state space. 

IV. NUMERICAL CALCULATIONS 

In this section we present results of numerical calcula- 
tions that demonstrate the properties of the system in the 
different parameter regimes. In particular, we perform 
calculations for the resonant case, the high-frequency- 
oscillator case and the high-frequency-qubit case. We 
also vary the qubit bias e, or alternatively the angle 9, in 
order to analyze its effect on the properties of the system. 

A. Energy-level spectrum 

In Fig. 13] we plot the energies of the lowest ten levels as 
a function of the coupling strength A in the resonant case 
huiQ = Eq. When e = A = 0, the ground state is non- 
degenerate and each higher energy level is doubly degen- 
erate. The separation between the levels is hujo, which 
is also equal to Eq. As A increases, the energy levels 
shift up or down, and several avoided crossings are en- 
countered. In the large A limit, all energy levels become 
doubly degenerate (i.e. they form pairs), including the 
ground state. The separation between the different pairs 
of energy levels in this limit is again huo- These results 
agree with the picture of the effective double- well poten- 
tial of Sec. nil CI For a small but finite bias e (i.e. small 
but finite 9) and small coupling strength A, the overall 
energy level structure is similar to that in the e = case, 
except that the levels do not approach each other as much 
at the avoided crossings. In the large A limit, there are 
no degeneracies: the energy levels are separated by the 
alternate distances e and hujQ — e. This structure reflects 
the small asymmetry in the double-well potential caused 
by a small tilt. For large 9 (i.e. sin 6* ~ 1), all features 
in the spectrum are suppressed, except for the overall 
decrease in the energy with increasing A. 

In order to examine the strong-coupling limit more 
closely, in Fig. |4] we plot the energy- level separation be- 
tween the lowest two energy levels. The results agree with 
the predictions of Eq. ([T^ : deep in the strong-coupling 
regime, the separation within the pairs of energy levels 
is given by E2n+2 - £^2n+i ~ (A/^iwo)" exp{-2{X/TiWQ f} 
for e = and by e for e ^ 0. 

In Fig. [5] we plot the energies of the lowest ten levels 
as a function of A in the case of a high-frequency oscil- 
lator, i.e. when Eq ^ TiuJo- As explained in Sec. IIII 51 if 
one considers a pair of energy levels, e.g. the lowest two 
energy levels, one has a modified effective qubit Hamilto- 
nian. When A = 0, one recovers the bare qubit Hamilto- 
nian. As A increases, the effective qubit gap A decreases 
and approaches zero in the limit A/ (S-Wq) — ^ oo. In Fig.lS] 
we plot the separations within the four lowest pairs of 
energy levels. The effective gap follows the shape of a 
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FIG. 3: (Color online) Lowest ten energy levels in the resonant 
case, i.e. when hujo/Eq — 1. The rescaled energy En/{Tiuio) 
with n = 1, 2, 10 is plotted as a function of the rescaled cou- 
pling strength X/{huuo)- Panels (a), (b) and (c) correspond to 
— 0,n/6 and 7r/3, respectively [recall that 6 — arctan(e/A)]. 

Gaussian function times a Laguerre polynomial, vanish- 
ing at the zeros of the Laguerre polynomial. As 9 is 
increased from zero, i.e. as the ratio A/e decreases, the 
dependence of the energy-level separation on the coupling 
strength becomes weaker (This phenomenon can be seen 
by comparing the different panels in Fig. Note that 
the location of the peaks does not change, but the ef- 
fect of the gap on the energy levels becomes smaller with 
increasing 9. 

In Fig. [7] we plot the energies of the lowest ten levels 
as a function of A in the case of a high-frequency qubit, 
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FIG. 4: (Color online) The separation between the lowest two 
energy levels in the resonant case, i.e. when huo/Eq = 1. In 
(a), the rescaled energy separation (E2 — E\)/{fiu)o) is plot- 
ted as a function of the rescaled coupling strength \/(hu)({). 
The blue, solid line corresponds to S = 0; the green, dashed 
line corresponds to ^ = 7r/6; and the red, dash-dotted line 
corresponds to = 7r/3. In (b), the same data is plotted 
on a logarithmic scale in order to make a comparison with 
the formula E2 - Ei = A exp{-2(A/;ia;o)^} from Eq. (fTSll : 
X and y in the axis labels refer to the axis labels in (a). 
The black, dotted line shows the asymptotic behavior of the 
above formula. The good fit between the blue and black lines 
means that the numerical results agree with the results of 
Sec, mil (We could extend the range of agreement by plotting 
log[log(A/?icjo) — log(j/)]; however, we are mostly interested in 
demonstrating the agreement for large values of x, where this 
modification would have little effect on the shape of the blue 
curve). Similar figures can be generated for the other values 
of hujo/Eq. However, we do not show such figures here. 



i.e. when Eq ^ Hloq. The most dramatic effects occur 
for 9 = 0. The ground-state energy remains essentially 
constant between A = and A = -/H(3oS/2. Beyond this 
point the ground-state energy decreases indefinitely with 
increasing A. Furthermore, below the critical point, the 
low-lying energy levels approach each other with increas- 
ing A as if they were going to collapse to one point, as 
would be expected for a vanishing lLiq. Above the critical 
point, the energy levels form pairs whose intra-pair sep- 
aration decreases with increasing A. The above scenario 
is suppressed as 9 is increased. There is no longer any 
sign of a critical point, and the energy-level separations 
are independent of A. 
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FIG. 5: (Color online) Lowest ten energy levels in the case of 
a high-frequency oscillator; fiujo /Eq = 10. Here we only show 
the results for 6 = 0, because the overall appearance of the 
plots is independent of 9. More details can be seen in Fig. [S] 



B. Squeezing, entanglement, and 'cat-ness' in the 
ground state 

One obvious possibility for the preparation of squeezed, 
entangled or Schrddinger-cat states in the case of ultra- 
strong coupling is to have a ground state that exhibits 
one of these unusual properties. With this point in mind, 
in this section we analyze the oscillator's squeezing and 
cat-ness as well as the qubit-oscillator entanglement in 
the ground state for different choices of system parame- 
ters. 

As a first step, we plot the Q function and the Wigner 
function of the oscillator's state in the ground state of 
the coupled system. The Q function is given by 

Q{X,P) = -{X + iP\p,sc\X + iP), (31) 

TT 

where posc is the oscillator's reduced density matrix af- 
ter tracing out the qubit from the ground state, posc = 
Trq{|5'Gs) (^Gsll with I^^gs) being the ground state of 
the combined system, and the bra and ket in the above 
formula are coherent states: 



= exp{c 



a*d} I 



(32) 



The state |0) represents the vacuum state with the oscil- 
lator in its ground state. The Wigner function is given 
by 



Po 



X--X') e'P'^'dX' 



.(33) 

where the bra and ket are now eigenstates of the position 
operator x, i.e. they are highly localized in configuration 
space. The Q and Wigner functions for a sequence of A 
values are shown in Fig. |8l 

Going beyond the pictorial description shown above, 
a quantifier for both the squeezing and cat-ness of the 
oscillator's state is the set of two squeezing parameters 
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FIG. 6: (Color online) The rescaled energy separation 
{E2n+2 — E^n+i) / Eq witMn the lowest four pairs of energy 
levels [i.e. for n — (blue, solid line), 1 (green, dashed line), 2 
(red, dash-dotted line) and 3 (black, dotted line)], as a func- 
tion of A/(?icJo)- As in Fig.O we take fuoo/Eq = 10. In panels 
(a), (b) and (c), = 0, 7r/6 and 7r/3, respectively. Note that 
the minimum value on the j/-axis is given by sin0 and is dif- 
ferent in the three panels. 



in the x and p quadratures as well as the product of 
the quadrature variances (Note here that the oscillator's 
state is always mirror-symmetric with respect to the x 
axis in the setup under consideration, giving (p) = 0). 
After the appropriate conversion into dimensionless vari- 
ables, these quantifiers are given by 



(( 



= 4 X-(X] 
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FIG. 7: (Color online) Lowest ten energy levels in the case 
of a high-frequency qubit; huJo/Eq — 0.01. In panels (a), 
(b) and (c), 9 = 0, 7r/6 and 7r/3, respectively. The insets 
show enlarged views of the weak-coupling and strong-coupling 
regions. 



P-(P) 



K 



({x^{x)f)({p~{p)y 

— (l-f S,)(l-f Sp). 



(34) 



- 1 



The parameter K is equal to Ji^/A for a minimum- 
uncertainty state (including both coherent and 
quadrature-squeezed states) and is larger than that 
lower bound for any other state (including Schrodinger- 
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FIG. 8: (Color online) The Q function (left) and the Wigner 
function (right) of the oscillator's state in the ground state of 
the combined system. Here we take hujo/IS. = 0.1 and e = 0. 
The different panels correspond to A/(?!,a;o) = 0.5 (a,b), 2 
(c,d), 2.5 (e,f) and 3.5 (g,h). For clarity, we adjust the color 
scheme in the different panels such that the highest point is 
always black. The red and yellow colors also correspond to 
positive values. The white color corresponds to zero value. 
The blue color represents negative values of the Wigner func- 
tion. The oscillator's state goes from a coherent state with 
no photons (i.e. the vacuum state) in the absence of coupling, 
to a squeezed state for low to moderate coupling strengths 
and then to a qubit-oscillator entangled state for very strong 
coupling. Note that the state in panels g and h is highly non- 
classical, in particular highly entangled, even though this fact 
cannot be seen in the Q and Wigner functions. 



cat and qubit-oscillator entangled states). In Fig. [9] we 
plot the momentum squeezing parameter as a function 



FIG. 9: (Color online) The momentum squeezing parameter 
Sp as a function of X/{tnJo) for hujo/A — 0.1 (a), 1 (b) and 10 
(c). The different curves correspond to e/A = (blue, solid 
line), 0.1 (green, dashed line), 0.5 (red, dash-dotted line) and 
1 (black, dotted line). The oscillator's state becomes squeezed 
as the coupling strength A increases, but then it reaches a 
maximum and goes back to zero as the qubit and oscilla- 
tor get entangled in the strong-coupling regime. Note that 
the maximum achievable squeezing decreases with increasing 
hujo/A. 



of the coupling strength A. For small values of A, 
the squeezing increases with increasing A. However, 
as A increases further and the ground state becomes 
more and more entangled, the squeezing is lost. The 
maximum achievable squeezing is largest for the case of 
a high-frequency qubit, huo <C A. Indeed, as explained 
in Sec. IIII CI the oscillator's effective potential becomes 
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flatter and flatter as one approaches the critical point, 
leading to a momentum squeezing parameter close to 
-1. 

As \sp\ increases, one can ask whether the oscillator's 
state is a quadrature-squeezed, minimum-uncertainty 
state or it deviates from this ideal squeezed state. The 
answer to this question can be obtained by analyzing the 
parameter K. We do not show any plots of this param- 
eter here. The main results are as follows: For the case 
e = 0, K increases slowly and remains close to h^/A as 
Sp increases, but near the maximum squeezing point, K 
starts increasing rapidly and diverges for X/{hujQ) — > oo. 
For finite values of e, if increases slightly above ?i^/4, 
but then turns and goes back to as Sp goes back to 
zero in the strong-coupling limit. 

We have seen that squeezed states are obtained for 
weak to moderate coupling. The question now is what 
states we have for strong coupling. The Q functions 
and the Sp and K results discussed above do not dis- 
tinguish between a Schrodinger-cat state in the oscilla- 
tor and a qubit-oscillator entangled state. The Wigner 
function has negative values for moderately strong cou- 
pling (Figs. IHJl andHf), indicating nonclassical states of 
the Schrodinger-cat type (Note that quadrature-squeezed 
states have nonnegative Wigner functions). In order to 
distinguish more clearly between Schrodinger-cat states 
in the oscillator and qubit-oscillator entangled states, we 
now analyze the entanglement properties in the ground 
state. 

The entanglement is quantified by the entropy S of the 
qubit's state. This quantity is obtained by calculating the 
ground state of the combined system I^'gs), using it to 
obtain the qubit's reduced density matrix in the ground 
state pq = Trosc{|*Gs) (^'gs|}, and then evaluating the 
entropy of that state S = — Tr{pqlog2 Pq}. 

In Fig. [To] we plot the qubit's ground-state entropy as a 
function of A. For e = the entropy increases from zero 
to one as A increases from zero to values much larger 
than all other parameters in the problem. Demonstrat- 
ing the fragility of the entangled states in the large- A 
limit. Fig. [TU] shows that the entanglement drops rapidly 
(especially for large values of A) when e is increased. 

By comparing Figs. 1^ and [TUl we can see that the rise 
in the qubit-oscillator entanglement is correlated with 
the reversal of the squeezing. One therefore goes from a 
squeezed state in the oscillator to a qubit-oscillator en- 
tangled state. We do not find any set of parameters where 
the ground state contains an unentangled Schrodinger- 
cat state in the oscillator. 

The numerical results show that the case huo <C Eq is 
most suited for the preparation of squeezed states, as can 
be seen by comparing the maximum achievable squeez- 
ing in the different parameter regimes. The opposite case 
{hujo ^ Eq) is most suited for the preparation of entan- 
gled states, as seen from the extreme fragility of these 
states for the case TiuJo <^ Eq. In fact, all the nonclas- 
sical properties of the ground state are suppressed as e 
is increased from zero to values larger than A. We shall 
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FIG. 10: (Color online) The qubit's entropy S (which quan- 
tifies the qubit-oscillator entanglement) in the ground state 
as a function of \/{hujo). The ratio fiujo/A is 0.1 in (a), 1 in 
(b) and 10 in (c) , and the different curves in each panel corre- 
spond to e/A = (blue, solid line), 0.1 (green, dashed line), 
0.5 (red, dash-dotted line) and 1 (black, dotted line). For 
e = the qubit-oscillator entanglement increases from zero to 
one as A is increased, regardless of the relation between Tiloo 
and Sq. However, the entanglement drops rapidly (especially 
for large values of A) as e is increased, i.e. when the qubit is 
moved away from the degeneracy point. 



come back to this point in Sec. I VII 

In Fig. [11] we examine the value of A at which the 
qubit's ground-state entropy has the values 0.1 and 0.5. 
These curves serve as indicators for the onset of qubit- 
oscillator entanglement, which is related to the instabil- 
ity encountered in the semi-classical calculation. For a 
high-frequency qubit, the sharp rise in entanglement oc- 
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FIG. 11: (Color online) The value of A/A at which the qubit's 
ground-state entropy has the values 0.1 (red, lower solid line) 
and 0.5 (blue, upper solid line) plotted as a function of hujo/A 
(Note the logarithmic scale). Here we take e = 0. The straight 
lines are given by the formulae A — y^hujoEq/2 (dotted line) 
and A — Tiujo (dashed line), which we have obtained in Sec. lIIII 
For small values of hujo/A, the onset of entanglement occurs 
when A = For large values of ftcJo/A, the onset 

of entanglement occurs when A ~ hujo, in agreement with the 
dependence explained in Sec. IIIII 



curs at A = y''hujQEq/2 which agrees with the instabihty 
condition of Sees. IIII CI and IIIIDI For a high-frequency 
osciUator, the entanglement rises to large values when 
A ^ hcdQ, in agreement with the analysis of Sec. IIII Bl 



V. PREPARATION AND DETECTION OF 
NONCLASSICAL STATES THROUGH IN-SITU 
PARAMETER AND STATE MANIPULATION 

We have seen in Sec. llVl that oscillator squeezed states 
and qubit-oscillator entangled states can occur natu- 
rally as ground states of the strongly coupled system. 
Schrodinger-cat states in the oscillator, i.e. not involv- 
ing entanglement with the qubit, do not occur as ground 
states of this system. 

One method that has been proposed for the generation 
of oscillator Schrodinger-cat states in the context of cav- 
ity QED [2] can be considered here as well: One prepares 
a qubit-oscillator entangled state of the form 

l={\a)(^\q,) + \~a)(g>\q2)), (35) 

with and \q2) being any two orthogonal qubit states 
and |±a) being coherent states of the oscillator with a 
large value of |a|. One now measures the qubit in the 
(|gi) ± \q2))/V^ basis. Depending on the outcome of the 
measurement, the state of the oscillator is projected into 
one of the two states 

4(l«)±|-«)), (36) 



each of which is a Schrodinger-cat state. Since the ground 
state well above the critical point is approximately given 

by 

-^(|a)®|i) + |-a)®|t)), (37) 

with a = xq, the above procedure could also be imple- 
mented in the system under consideration (We shall also 
give an alternative procedure below). Hence all three 
types of nonclassical states that we consider in this pa- 
per can be generated in principle. 

One important question that arises in the case of strong 
qubit-oscillator coupling is whether it is possible to de- 
tect the nonclassical states in spite of the always-present 
strong coupling. The answer is yes, in principle. An im- 
portant point to note in this context is that, as shown 
in Sees. Illll and IIVI the energy eigenstates of the system 
are approximately product states when the qubit is bi- 
ased far from the degeneracy point, i.e. for large values of 
e. One could therefore say that, for certain procedures, 
the qubit and oscillator can be made to effectively decou- 
ple from each other by biasing the qubit away from the 
degeneracy point. 

We now discuss some possible procedures for the exper- 
imental observation of the different nonclassical states. 
Since all three types of nonclassical states occur in the 
case of a high-frequency, adiabatically adjusting qubit, 
we focus on this case. An experiment could start by 
preparing the ground state with the qubit biased at the 
degeneracy point. Under suitable conditions, this step 
could be achieved by biasing the qubit at the degeneracy 
point and letting the system relax to its ground state. 
The state would then be either a squeezed state or an en- 
tangled state, depending on the coupling strength. One 
could then move the qubit away from the degeneracy 
point for measurement purposes. If the change is slow on 
the timescales of both the qubit and the oscillator, the 
system will remain close to its ground state, adiabatically 
following the bias-point shift (There might be substan- 
tial excitation out of the ground state in the double- well 
regime where the separation between the lowest energy 
levels is small; however, this situation will come to an end 
when the double-well potential transforms into a single- 
well potential). If the system follows its instantaneous 
ground state during the bias-point sweep, the nonclas- 
sical state will be lost. The sweep therefore has to be 
fast at least compared to the period of the oscillator. If 
the sweep is adiabatic with respect to the qubit but fast 
with respect to the oscillator, the qubit will follow the 
change adiabatically, while the oscillator will be frozen 
in its initial state. In the case where the initial ground 
state is a squeezed state, one would achieve the effective 
qubit-oscillator decoupling while preserving the squeezed 
state for the measurement step of the experiment. In the 
case where the initial ground state is (approximately) the 
entangled state given in Eq. p7l) . the qubit will end up 
in its ground state at the final bias point, and this state 
will be independent of the state of the oscillator. As a 
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result, the oscillator is left in a Schrodinger-cat state. If 
the bias-point sweep is fast on both the qubit and oscil- 
lator timescales, both subsystems will be frozen in their 
initial state during the sweep, such that one ends up with 
an entangled state at the end of the sweep. 

The state of the oscillator can be reconstructed using 
Wigner tomography, which could be implemented follow- 
ing the experiment in Ref. [2l| . In that experiment, the 
oscillator was put into resonance with a qubit that was 
initialized in its ground state, and the excitation proba- 
bility of the qubit as a function of time was determined by 
performing an ensemble of measurements. By decompos- 
ing the signal into its Fourier components, it was possi- 
ble to extract the occupation probabilities of the different 
photon- number states. When combined with the ability 
to shift the oscillator's state (through the application of 
a classical driving signal) before the measurement, full 
Wigner tomography becomes possible. In the case of a 
low-frequency oscillator, the transfer of excitations be- 
tween the qubit and the oscillator can be induced by 
driving the red or blue sideband, as was done in the ex- 
periment of Ref . ^3,]. Since, as explained above, the qubit 
is effectively decoupled from the resonator when e >• A, 
the exchange of excitations between the qubit and the 
oscillator would not be efficient at the measurement bias 
point, which seems to pose a dilemma for the proposed 
experiment. This difficulty can be circumvented, how- 
ever, by using a second, weakly coupled qubit for mea- 
surement purposes. 

Wigner tomography can be used to demonstrate 
squeezed and Schrodinger-cat states in the oscillator. 
Qubit-oscillator entangled states could be demonstrated 
by measuring the correlation between the states of the 
qubit and oscillator. Starting from the ground state, if 
the qubit is measured and found to be in the state |t), the 
oscillator must be in the state |— a), with a = xq. If the 
qubit is found to be in the state |4-), the oscillator must 
be in the state \a). The observation of only these cor- 
relations, however, is not sufficient in order to establish 
the presence of quantum correlations. For that purpose 
one has to perform measurements in more than one set of 
bases. The additional qubit basis can be (|t) ± \i))/V^: 
If the qubit is found to be in the state + \\.))/^/2, the 
oscillator must be in the state (|q;) -I- \—a))/^/2, and a 
similar relation holds for the minus signs. The two states 
± \—a))/^/2 can be distinguished through the fact 
that the state with the plus sign contains only even pho- 
ton numbers while the state with the minus sign contains 
only odd photon numbers. 

Finally it should be noted that after the bias-point 
sweep, the resulting state would not be a stationary state 
and would therefore have a time dependence. This time 
dependence has to be taken into account in the measure- 
ment sequence. Furthermore, when the qubit is biased 
such that it is in one of the two ct^ eigenstates, say |t), 
the effective oscillator potential will be shifted from the 
point X = 0, and one must take into account this shift 
when analyzing the post-sweep dynamics. 



VI. DECOHERENCE 

We now turn to the question of how coupling to the en- 
vironment affects the prospects of preparing and observ- 
ing nonclassical states in the system under consideration. 
Following a standard procedure 28, 32, 33], we analyze 
the effects of the environment by first determining the 
energy eigenstates of the system in isolation and then 
analyzing the relaxation and dephasing rates that gov- 
ern the decoherence between the different energy eigen- 
states. For our purposes it will be sufficient to consider 
Markovian decoherence dynamics. 

Since the preparation of nonclassical states above re- 
quired biasing the system at the point e = 0, we focus 
on this case. In Sees. IIIII and IIVI we found two types 
of low-lying energy eigenstates, depending on the cou- 
pling strength. For small values of A, where the ground 
state involves a squeezed state of the oscillator, the en- 
ergy eigenstates are slightly modified from the energy 
eigenstates in the absence of qubit-oscillator coupling. 
For large values of A, the low-lying energy eigenstates 
are superpositions similar to that given in Eq. (l37l) . As 
we shall see below, there is a qualitative difference in how 
these two types of states are affected by the environment. 

The relaxation rate Ti^j and the dephasing rate T^^ij 
involving states i and j are given by [2a . [s^ 



r^,,, = nSiO) X mA\i)~{j\A\j) 



(38) 



where S{uj) is the environment's spectral density of the 
relevant environment operator at frequency uj, and A is 
the system operator through which the system couples 
to the environment. 

In order to go further with the analysis, we need to 
specify the operator A that describes the coupling to the 
environment; there is one such operator for each decoher- 
ence channel. It was mentioned in the introduction that 
the availability of the tuning parameter e can be seen 
as an advantage of solid-state qubits in comparison to 
natural atoms in cavity QED. On the other hand, how- 
ever, this property also opens an additional channel for 
noise and the environment to couple to the system. The 
operator associated with the parameter e is 172, and the 
coupling through this operator is typically the most detri- 
mental for superconducting qubit circuits. We therefore 
start by considering this decoherence channel. 

For coupling through the operator ci^, the relaxation 
and dephasing rates are proportional to the quantities 
and \{i\az\i) - {j\ respectively. For 

e = and small A (and avoiding the resonant case), the 
energy eigenstates are approximately product states with 
the qubit being in an eigenstate of (t,j, , to which we refer 
as |±), and the oscillator having a certain number of 
photons n. The above quantities are then approximately 
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given by the corresponding values for the qubit states: 

(n,+|(73 |n',+) = |n',-) =0. (39) 

These expressions suggest that the system relaxes with 
a rate equal to that of the isolated qubit. The vanish- 
ing of the dephasing rate in this approximation has the 
same origin as its vanishing for an isolated qubit at the 
degeneracy point, namely the fact that the energies are 
insensitive to variations in e to first order. This prop- 
erty points out an important requirement for the above 
expressions to be valid: the energy-level separation must 
be much larger than the transverse fluctuations in the 
Hamiltonian. These fluctuations, which are transverse 
at the degeneracy point, are responsible for dephasing 
away from the degeneracy point. As a result, in order 
to obtain the degeneracy-point protection from dephas- 
ing, the energy-level separation must be large compared 
to the dephasing rate when the latter is calculated away 
from the degeneracy point: 



r^,^j = ^5(0) X |(n,t|<T. |n,t) 



(n',i|<7,|n',;)r 



4775(0). 



(40) 



This parameter is typically of the order of 100 MHz, 
which corresponds to a dephasing time of 10 ns. At the 
degeneracy point, an isolated qubit is protected from this 
noise because the qubit 's gap is typically larger than 1 
GHz, and the qubit-environment coupling is transverse to 
the qubit Hamiltonian. Similarly, mildly squeezed states 
(whose energy-level separations are comparable to those 
of the simple product states of the uncoupled system) 
should be protected from dephasing caused by the weak, 
transverse coupling to the environment. The situation is 
quite different for large values of A, where the low-lying 
states are highly entangled states with a very small sep- 
aration within an energy-level pair. For example, when 
e = 0, the lowest two energy eigenstates are given by 



1 



V2 



{\a)®\\)±\-a)®m. 



(41) 



with a = xq. The energy separation between these states 
can be obtained from the WKB calculation of Sec. IIII CI 
When e exceeds this (possibly very small) energy sep- 
aration, the energy eigenstates are simply the product 
states \a) ® ||) and |— a) ® W). In order for the entan- 
gled states to be robust against fluctuations in e, their 
energy-level separation (at e = 0) must be much larger 
than the qubit's dephasing rate of about 100 MHz. Since 
the qubit and oscillator frequencies can be of the order of 
1 GHz, one could obtain an entangled ground state that 
is separated from the first excited state by 100 MHz or 
more. For example, taking A = Tiloq = A = 1 GHz, we 
obtain a qubit ground-state entropy of 0.85 and an en- 
ergy separation of 138 MHz. However, it would be highly 
desirable to reduce the qubit decoherence rates, and in 
principle this should be possible in the future using better 
materials and circuit designs. 



Even though superconducting harmonic oscillators 
generally have much higher quality factors than super- 
conducting qubits, it is instructive to briefly discuss the 
effect of oscillator decoherence. The oscillator typically 
couples to its environment through the same operator 
that describes its coupling to the qubit. In the present 
problem, this operator is x, or equivalently a -t- a^. In a 
free oscillator, the relevant matrix element for purposes 
of determining the decoherence rates is given by 

(n| (a + a)) \n') \/n5n-i,n' + V"- + l(5„+i^„/. (42) 

Using the relations for the relaxation and dephasing rates 
in Eq. (p8)) . one finds that at low temperatures the effect 
of the environment is to cause decay to the ground state 
through the loss of photons; photons are lost one by one, 
with a rate that is proportional to the photon number 
(Note that since {n\ {a -f a^) \n) — 0, no pure dephasing 
occurs in this system). The photon- loss process is de- 
scribed by the jump operator a. In the strongly coupled 
qubit-oscillator system with the double-well effective po- 
tential, the effect of the environment will be drastically 
different from the simple photon-loss dynamics. Since 
the effective potential near each one of the local minima 
has the same shape as the free-oscillator potential, the 
energy eigenstates |n, ±) will be entangled states where 
the qubit is either in the state |t) or jj,), correlated with 
an oscillator wave function that is given simply by the 
free-oscillator state \n) shifted to the left or right by a 
distance xq. One therefore finds the matrix elements 



{n,±\{a + a))\n' ,±) = ^/n5n-i,n' + \/ n + l(5„+i^„' 
{n,+\{a + a^)\n',-) = 0. (43) 

The fact that the matrix elements in the first line co- 
incide with those of the free oscillator implies that the 
relaxation rate will be equal to the free-oscillator relax- 
ation rate. Note that this process is no longer described 
by the jump operator a, but rather by a properly shifted 
annihilation operator, a± y^muJo/{2h)xo, depending on 
whether one is dealing with the left or right well in the 
effective double-well potential. The relaxation process 
does not change the state of the qubit, or any superposi- 
tion involving the left and right wells. Even though the 
relation 



(n,±\ la- 



\n,±) 







(44) 



would suggest that no dephasing should occur between 
the different energy eigenstates, the alternative basis with 
localized states {|n,t) , gives 

I (n, tl (a + a^) \n, t) - (n, ;| (a + a^) \n, i) |' = ^^xl 

(45) 

This result implies that the coupling to the environ- 
ment will cause dephasing in any quantum superposi- 
tion involving the two wells, with a tendency to local- 
ize the wave function in one of the wells. The rate 
of this process will be proportional to the product of 
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the environment's power spectrum at zero frequency and 
the combination of matrix elements given above. This 
last quantity is proportional to Xq, and it grows indefi- 
nitely with increasing coupling strength. If we assume 
that the power spectrum at zero frequency is compa- 
rable to that at the oscillator frequency, the dephas- 
ing rate can be much larger than the decay rate of the 
free oscillator because of the largeness of the quantity 
|(n,t| (a + a'') |n,t) - {n,i\ (a + at) For exam- 

ple taking the oscillator's decay rate to be 0.1-1 MHz and 
considering states where a and are of typical size -v/S 
(i.e. about five virtual photons in the ground state), we 
find that the dephasing rate can be of the order of 10-100 
MHz, which is comparable to the (T^-mediated dephasing 
rate. 

We finally discuss the question of temperature. We 
have considered the possibility of preparing nonclassical 
states by letting the system cool down to its ground state. 
One must therefore make sure that the energy-level sep- 
aration between the ground state and the first excited 
state is larger than the ambient temperature. In super- 
conducting circuits, the temperature is typically around 
20 mK, which can be converted to roughly 1 GHz in 
frequency units. The ground state must be separated 
from the excited states by at least that amount in or- 
der to achieve high-fidelity preparation of the ground 
state. Typical qubit and oscillator frequencies are in 
the few-gigahertz range, not much higher than typical 
temperatures. Squeezed ground states in the oscillator 
should be separated from the excited states by an en- 
ergy comparable to the one present in the uncoupled sys- 
tem, implying that the preparation of these states should 
be possible. The entangled ground states that occur for 
strong coupling, however, are separated from the first ex- 
cited states by energy gaps that decrease rapidly with in- 
creasing qubit-oscillator coupling strength. If this energy 
gap becomes smaller than the 1 GHz temperature level, 
one would not be able to prepare the entangled ground 
state simply by letting the system cool down to such a 
state. However, one could let the system cool down to its 
ground state away from the degeneracy point and then 
adiabatically shift the bias point to one with an entangled 
ground state. Provided that the thermalization rate is 
sufficiently low, it is not necessary to have a degeneracy- 
point energy-level separation that is larger than the tem- 
perature. The 1 GHz temperature level should therefore 
not be seen as a fundamental obstacle to the preparation 
of entangled ground states. 



VII. CONCLUSION 

We have analyzed the properties of a strongly cou- 
pled qubit-oscillator system, focusing on the potential 
of this system for the preparation of nonclassical states. 
These states include squeezed states and Schrodinger-cat 
states of the oscillator, as well as qubit-oscillator entan- 
gled states. 



We have compared the predictions of four differ- 
ent analytical approaches: the weak-coupling approx- 
imation, the adiabatically-adjusting-oscillator approxi- 
mation, the adiabatically-adjusting-qubit approximation 
and the semiclassical calculation. Each one of these four 
approaches is well suited for analyzing the behaviour of 
the system under a certain set of assumptions. Thus 
the combination of the results provides a rather thor- 
ough understanding of the qubit-oscillator system in the 
regime of ultrastrong coupling. We have also presented 
results of numerical calculations that reinforce the results 
of the analytical derivations. These results demonstrate 
the nonclassical properties of the energy eigenstates, and 
especially the ground state, of the system. 

We have discussed various possible experimental pro- 
cedures for the preparation and observation of nonclas- 
sical states. All three types of nonclassical states that 
we discuss in this paper can be prepared and detected in 
principle. 

We have also analyzed the effect of coupling to the 
environment on the system. We have shown that the de- 
coherence dynamics of the coupled qubit-oscillator sys- 
tem can be qualitatively different from the decoherence 
dynamics of the qubit or oscillator in isolation. We 
have shown that nonclassical states, particularly highly- 
entangled states, are highly susceptible to changes or 
fluctuations in the bias parameters. These results lead 
to the conclusion that high degrees of control and low 
noise levels will be required for the preparation of robust 
nonclassical states in the ultrastrong-coupling regime. 
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tional Science Foundation (NSF) grant No. 0726909. 

Note added in proof: Recently ultrastrong coupling 
between a superconducting fiux qubit and a coplanar 
waveguide resonator has been demonstrated [33 |. 

Appendix A: Oscillator's effective kinetic energy 
in the case of a high-frequency qubit 

In this appendix we briefiy discuss the modification 
to the oscillator's kinetic energy in the case of a high- 
frequency qubit (see Sec. IIII C"|) . One could perform the 
calculation by considering only the kinetic-energy term, 
without considering the specific form of the trapping po- 
tential, i.e. by using the completely delocalized momen- 
tum eigenstates as a basis for the calculation. Such a 
calculation, however, turns out not to lead to simple, 
transparent results. We therefore consider the correc- 
tions that one would need to add to the kinetic-energy 
term in the effective oscillator Hamiltonian starting from 
the eigenstates of the free-oscillator Hamiltonian. 

The eigenstates of the Hamiltonian iJho can be ex- 
pressed in the position basis as 

|fc) = / dxMx)\x). (46) 
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Taking into account the high-frequency, adiabaticaUy ad- 
justing qubit, and for definiteness taking the case of the 
qubit's ground state, the above eigenstates are modified 
as foUows: 



k) = I dxipk{x) \x) <Si \g{x)) 

dxMx)Y.^*k'i^)\k')(^\9ix)), (47) 



where \g{x)) has the same meaning as \gx) from Sec. lIII Cl 
One can now obtain the correction to the kinetic-energy 
term as: 



P 

2m 



^ 10 = y j dxidx2 ^l{xi)^k' jxi] 

k'V 



rv{x2)Mx2) (fc'l ^ 10 ((.9(2:1)1 g{x2)) 1) .(48) 
The factor between parentheses represents the relative 



correction to the kinetic-energy term in the Hamihonian. 
This factor can be estimated as: 

(5(.i)|,(x.))-l^cos(^ii^l^)-l 

[ip{xi) - lf{x2)f 



g{xi ~ X2) 



muJoE^ ■ 



(49) 



In deriving this expression we have taken the case e/A <C 
1 (for which the relative correction is maximum) and 
taken x to be of the order of the characteristic oscilla- 
tor length. By comparing the above expression to the 
relative correction in the potential-energy term found in 
Sec. lIII Cl i.e. g"^ /{mujQEq), one can see that the kinetic- 
term correction is negligible when huio <C Eq. 
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